Dissipative dynamics of a charged particle in the field of three plane waves: chaos and 

control 
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The chaotic dissipative dynamics of a charged particle in the field of three plane waves is theoret- 
ically (Melnikov's method) and numerically (Lyapunov exponents) investigated. In particular, the 
effectiveness of one of such waves in controlling the chaotic dynamics induced by the remaining two 
waves is theoretically predicted and numerically confirmed. Two mechanisms underlying the chaos- 
suppression scenario are identified. One mechanism requires chaos-inducing and chaos-suppressing 
waves to have both commensurate wavelengths and commensurate relative (with respect to the 
remaining third wave) phase velocities, while the other mechanism allows the chaotic dynamics to 
be tamed when such quantitites are incommensurate. The present findings may be directly applied 
to several important problems in plasma physics, including that of the chaos-induced destruction of 
magnetic surfaces in tokamaks. 
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The interaction of charged particles with a wave packet 
is a basic and challenging problem appearing in many 
fundamental fields such as astrophysics, plasma physics, 
and condensed matter physics [1,2], to name just a few. 
As is well known, the structure of the wave packet reflects 
the particular physical situation that is being considered 
in each case. In this regard, a classical, widely studied 
particular case is that of an infinite set of electrostatic 
waves having the same amplitudes and wave numbers, 
zero initial phases, and integer frequencies [3]. While 
the Hamiltonian approach to the general problem is suit- 
able in many physical contexts [4], the consideration of 
dissipative forces seems appropriate so as to take into ac- 
count such diverse phenomena as electron-ion collisions, 
energy losses through synchrotron radiation, and stochas- 
tic heating of particles in plasma physics. In any case, 
stochastic (chaotic) dynamics already appears (can ap- 
pear) when the wave packet solely consists of two plane 
waves [5-8]. As is well known, this non-regular behavior 
of the charged particles may yield undesirable effects in a 
number of technological devices such as the destruction 
of magnetic surfaces in tokamaks [9]. Thus, in the con- 
text of plasma physics inter al. , it is natural to consider 
the problem of regularization of the chaotic dissipative 
dynamics of a charged particle in a wave packet by a 
small amplitude uncorrelated wave, which is added to 
the initial wave packet. This Letter studies the simplest 
model equation to examine this problem: 

X -\- ^x — [Eq sin ikox — Lo^t) + Ec sin ikcX — ujct)] 
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where the amplitudes Eq, Ec, Es, wave numbers ko,kc,ks, 
and frequencies ujq,lUc,uJs correspond to the "main", 
chaos-inducing, and chaos-suppressing waves, respec- 
tively, '^ is an initial phase, e and m are the charge 
and mass of the particle, respectively, and where 
weak dissipation (7 <ti 1) and non-uniform amplitudes 



{Ec,s/ Eq < 1) are assumed. In a reference frame moving 
along with the main wave, Eq. (1) transforms into the 
equation 
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where ^ = kox — cjoi, f^o = {skoEo/m) ' , t 
6 = 7Wo^o^' V = 7^0 \ £c,s = Ec,s/Eo, Kc,s = 
kc,s/ko, and ilc^s = {^c,s ~ ^okc.s/ko)n.Q are all di- 
mensionless variables and parameters. Also, flg/Qc = 
KsVos/ (KcVqc) = XcVos/iXsVoc), where vqc {vqs) and 
Ac (As) are the relative phase velocity with respect to 
the main wave {voc,s = ^c,slkc,a — too/ko) and the wave- 
length of the chaos-inducing (-suppressing) wave, respec- 
tively. The parameters kQ,u!o, Eq, and fio are held con- 
stant throughout. Physically, Eq. (2) can be used to de- 
scribe a number of important problems in plasma physics: 
the dissipative dynamics of a charged particle in the field 
of three electrostatic plane waves, the bounce dissipa- 
tive motion of a charged particle trapped in a toroidal 
magnetic field which is perturbed by two (one chaos- 
inducing and the other chaos-suppressing) superposed 
electrostatic waves, and the suppressory effect of a (sec- 
ondary) resonant magnetic perturbation on the destruc- 
tion of magnetic surfaces by a (primary) resonant mag- 
netic perturbation in a tokamak. Since Eq. (2) represents 
a perturbed pendulum (0 < 6, 77, ec,Ss ^ 1), one can ap- 
ply Melnikov's method (MM) [10,11,2] to obtain analyt- 
ical estimates of the ranges of parameters (Eg, kg, ujg, ^) 
for suppression of the chaos existing in the absence of the 
chaos-suppressing wave. The application of MM to Eq. 
(2) gives [12] the Melnikov function 



M"^ (to) = -D"^ + A^ sin [VLcTq] + B^ sin {Q.sTq + *) , 

(3) 
with D± = 877 ± 21:5, A^ = Ae^L^ {Kc,Vlc), 
S± = AesL^{Ks,VLs), and L^{K,n) = 



± /g°° sech T cos [2K arctan (sinh r) =F ^t] dr (see Fig. 
1), where the positive (negative) sign refers to the top 
(bottom) homochnic orbit of the underlying conservative 
pendulum. It is straightforward to obtain the follow- 
ing properties: (i) L^ {K,±Q,) = —L^{K,^fl); (ii) 
L+ (K, O) vanishes completely in the region (K, 51 < 0) 
but in a very narrow neighborhood of the line 51 = 0; 
(iii) There exists a particular relationship, fl = aK with 
a ~ (1 + \/oj /2 (i.e., the golden ratio), which deter- 
mines the direction of the main tongue-like region where 
L^ (K,n) presents its largest values (i.e., for small 
K values); (iv) L^ {K ^ 0,n) = (7r/2) sech(7rn/2), 
i±(X = n,f7 = 0) = 0,n = 1,2,... ; and {v) 
L^{K,VL) e [-7r/2,7r/2], \I{K,VL). It is well known 
that the simple zeros of the Melnikov function imply 
transversal intersection of stable and unstable manifolds, 
giving rise to Smale horseshoes and hence hyperbolic 
invariant sets [11]. The mechanism for taming chaos 
by an added weak periodic perturbation is then the 
frustration of a homoclinic (or hctcroclinic) bifurcation 
[13]. Consider first the case with no chaos-suppressing 
wave {Es =0). With fixed 6,r],ec, one sees that for 
sufficiently small values of kc the best chance for the 
occurrence of a homoclinic bifurcation, i.e.. 



\2r]±Tr6/2\ 



< \L^{K„n, 



(4) 



takes place when \voc\ is close to |woc,max|, where 
I'Oc.max = ±a51o/fco is thc most chaotic relative phase 
velocity (cf. properties (i) and (iii)). Figure 2 shows an 
illustrative example of this prediction, where the MM- 
based approximation (top panel) is compared with Lya- 
punov exponent (LE) calculations of Eq. (2) (bottom 
panel) . One sees that the largest values of the maximal 
LE (black points) lie inside a narrow tongue-like region 
which is near the theoretical estimate n = aK (yellow 
line in Fig. 2). Note that the maximal (positive) LE 
increases as Kc is increased, which is coherent with thc 
corresponding growth of the chaotic threshold function 
|Z/+ {K,V,)\. With fixed S,ri,ec, chaotic motion is possi- 
ble when ric = (i.e., when the relative phase velocity 
vanishes) and kc ^ nko , n positive integer (cf. property 
(iv), although this possibility decreases as Ac -^ 0. Also, 
as expected, property (i) means that the chaotic thresh- 
old depends on thc absolute value of thc relative phase 
velocity, but not on its sign (cf. Eq. (4)). Let us sup- 
pose in thc following that, in the absence of any chaos- 
suppressing wave {Eg ~ 0), the associated Melnikov func- 
tion Mq (to) = —D^ + A^ sm{n,cTo) changes sign at 
some To, so that thc charged particle exhibits (at least 
transient) chaotic behavior. To study the taming effect 
of the chaos-suppressing wave in the most favorable situ- 
ation, the case of a subharmonic {fls = P^c^P a positive 
integer) resonance is analyzed below. For this case, one 
has 






which permits one to identify two different physical situ- 
ations (or mechanisms) for the added chaos-suppressing 
wave to tame chaotic charged particles: I. Chaos- 
inducing and chaos-suppressing waves having both com- 
mensurate wavelengths (A^/Ac = m/n, m, n positive 
integers) and commensurate relative phase velocities 
(vos/voc = m'/n',m',n' positive integers) such that p = 
m'n/(n'm). II. Chaos-inducing and chaos-suppressing 
waves having both incommensurate wavelengths and in- 
commensurate relative phase velocities such that the quo- 
tient (wos/^oc) / (As/Ac) satisfies Eq. (5). 

Let ris = p^c, P a positive integer, such that the 
relationships * = ^lopt = 7r[4m -I- 3 - p (4n + l)]/2, 
7r[4m + 5 - p{4n -\- l)]/2, 7r[4m + 3 - p(4n-l)]/2, 
7r[4?Ti + 5 — p (4n — l)]/2 are independently satisfied for 
some positive integers m and n for the parameter re- 
gions {K,n) where [L± (i^T^, 51c) > 0,L±(i^s,^s) > 0], 
[L^iKc,nc) > 0,L^iKs,ns) < 0], [L^iKc^c) < 

0, L± {K,,n,) > 0], and [L± (Kc, 51c) < 0, L± {K,, n,) < 

01 , respectively. Then the frustration of a homoclinic bi- 
furcation occurs (i.e., M^ (tq) always has the same sign) 
if and only if the conditions 
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£r^.„ = (i-|z?±M±|)i?±, 

^ — D± / 2 

i?^ = ec\L^{Kc,nc)/L^iKs,ns 



(6) 



are fulfilled for each region {K,n), respectively. (The 
proof will be given elsewhere [12].) 

Now one can make the following remarks. First, 
for the above {K, 51) regions, this result requires hav- 
ing (* = 7r,7r/2,0,37r/2), (* = 0, 37r/2, tt, 7r/2), (* = 
0,7r/2,7r,37r/2), (* = tt, 37r/2, 0, 7r/2), respectively, for 
p — Am — 3, 4?Ti — 2, 4to — 1, 4to (m — 1, 2, ...), respec- 
tively in each case. Second, the effectiveness of the chaos- 
suppressing wave decreases as the resonance order p is 
increased (cf. Eq. (6)). This is relevant in experimental 
realizations of the control where the main resonance case 
is then expected to be the most favorable to reliably tame 
the chaotic dynamics. Third, for the main resonance case 
(51s = ^c), one has a more accurate and complete esti- 
mate [14] of the regularization boundary in the '^ — Ss 
parameter plane: 



=Fcos*± Wcos2* 



1 - {D±/A±f 



i?±. 



for B^ > and A^ , respectively, and 



(7) 



zbcosv]/ ± Wcos2 vj/ 



1 - {D±/A±y 



R^ 



(5) 



(8) 
for B^ < and A^ , respectively. The two signs be- 
fore the square root apply to each of the sign super- 
scripts of ef, which, in its turn, is independent of the 
sign of D^ . Also, the area enclosed by the boundary 



functions is given by 4 (|Z)^/A^|) i?^. The relevance of 
the above theoretical results on strange chaotic attractor 
elimination is confirmed by means of LE calculations of 
Eq. (1) [15]. Figure 3 shows the results corresponding 
to an illustrative example of the mechanisms of type I 
(top) and type II for the main resonance case {p = 1). 
In the absence of the chaos-suppressing wave {Sg = 0), 
Eq. (1) presents a strange attractor with a maximal 
LE A+ {es = 0) = 0.073 ± 0.001 (bits/s). The diagrams 
in this figure were constructed by only plotting points 
on the grid when the respective LE was larger than 
(cyan points) or than A"*" [ss = 0) (magenta points), and 
with the black dashed-line contour denoting the theoret- 
ical boundary function (cf. Eq. (7)) which is symmetric 
with respect to the optimal suppressory value '^opt = tt. 
Regarding the mechanism of type I, one typically finds 
that complete regularization (A+ (ss > 0) ^ 0) mainly 
appears inside the island which symmetrically contains 
the theoretically predicted area, while regularization by 
a type II mechanism seems to be insensitive to the ini- 
tial phase ^. In this latter case, the lowest value e^^^^ of 
the theoretical boundary function roughly coincides with 
the regularization threshold value. The regularization re- 
gions in the '^ — Ss parameter plane under the two types 
of mechanisms present different features. For type I, one 
typically finds that inside the regularization area which 
contains the predicted area, the two non-null LEs are 
identical and constant (A+ = A^ = —7/2), as is shown 



in the instance of Fig. 4 (top panel). This symmetry 
property of the contraction of the phase space volume for 
the predicted regularization areas in the ^ — £« parame- 
ter plane is an inherent feature of the type I mechanism. 
This property does not hold for the type II mechanism, 
since in this case the regularization region is not com- 
pletely " clean" (there exist isolated points corresponding 
to chaotic behavior) and furthermore the distribution of 
the non-null (negative) LEs is not perfectly uniform, as 
can be appreciated in the instance of Fig. 4 (bottom 
panel) . 

In conclusion, the onset of chaotic dissipative dynamics 
of a charged particle in the field of two plane waves and 
the effectiveness of an added plane wave in taming that 
chaos have been theoretically predicted and numerically 
confirmed. Two suppressory mechanisms were identi- 
fied: One mechanism requires chaos-inducing and chaos- 
suppressing waves to have both commensurate wave- 
lengths and commensurate relative phase velocities, while 
the other allows chaos to be tamed when these quantities 
are incommensurate. A more detailed discussion of the 
regularization routes by the two mechanisms for the main 
and higher resonances {p > 1) as well as of the Hamilto- 
nian limiting case will be given elsewhere [12]. It should 
be stressed that the results discussed in this work can be 
directly applied to the problem of the chaos-induced de- 
struction of magnetic surfaces in tokamaks besides other 
important problems in plasma physics. 
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Figure Captions 

Figure 1 (color online) . Function _L+ {K, il) (see tlie 
text). 



panel, kg = (\/5 — l) fcc/2, uqs = (-y/S + l) uoc/2. Black 
dashed-line contour indicates the predicted boundary (cf. 
Eq. (7)). System parameters: ko = luq = CIq = l,Cls = 
n^, Ec ^ 0.7', kc = 1.22, LUc = 2.26, 7 = 0.1. 



Figure 2 (color online). Top panel: Contour plot of 
the threshold function |L+ {K, 0)| with a gray scale from 
white (1.0 contour) to black (0.0 contour). Bottom panel: 
Grid of 200 x 200 points in the Kc — i^c parameter plane 
where cyan, magenta, and black points indicate that the 
respective LE was larger than 0.0, 0.07, and 0.14, respec- 
tively. System parameters: S = rj = 0.1, ec = 0.7. 

Figure 3 (color online). Grids of 100 x 100 points in 
the ^P — £5 parameter plane for the mechanisms of type 
I (top panel, kg = k^, Vqs = ^Oc) ^nd type II (bottom 



Figure 4. Top panel: Maximal LE A+ vs the ini- 
tial phase ^ for six values of the suppressory amplitude 
Es = 0.2,0.4,0.6,0.8,1.0,1.2 corresponding to the type 
I mechanism. Bottom panel: Maximal LE A+ vs sup- 
pressory amplitude £« for six values of the initial phase 
xp = 7r/3, 27r/3, tt, 47r/3, 57r/3 corresponding to the type II 
mechanism. Dashed-line contour indicates the predicted 
boundary (cf. Eq. (7)). System parameters as in the 
caption to Fig. 3. 
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